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Abstract 

This is the first part of a study of the quark propagator and the vertex function of the 
vector current on the lattice in the Landau gauge and using both Wilson-clover and overlap 
actions. In order to be able to identify lattice artifacts and to reach large momenta we use a 
range of lattice spacings. The lattice artifacts turn out to be exceedingly large in this study. 
We present a new and very efficient method to eliminate the hypercubic (anisotropy) artifacts 
based on a systematic expansion on hypercubic invariants which are not 50(4) invariant. A 
simpler version of this method has been used in previous works. This method is shown to 
be significantly more efficient than the popular "democratic" methods. It can of course be 
applied to the lattice simulations of many other physical quantities. The analysis indicates a 
hierarchy in the size of hypercubic artifacts: overlap larger than clover and propagator larger 
than vertex function. This pleads for the combined study of propagators and vertex functions 
via Ward identities. 



PACS: 12.38.Gc (Lattice QCD calculations) 



1 Introduction 

The study of the quark propagator and vertex functions has been extensively pursued in the 
hterature starting in the 70's [I]. Lattice QCD has more recently treated this issue A 
systematic treatment varying the quark actions has been followed by the CSSM collaboration [3]. 
The scalar part of the quark propagator is related via Ward identities to the pseudoscalar vertex 
function. The role of the Goldstone boson pole in the latter has been thoroughly discussed [1]. 

Leaving aside the latter issue, we will mainly concentrate on the vector part of the quark 
propagator, the one which is proportional to j). One of our main goals is to check the effect of the 

condensate which has been discovered via power corrections at large momenta to the gluon 
propagator and three point Green functions [3] -[HI- The values plotted in literature for Z^(p^) are 
extremely flat above 2 GeV At first sight this is a satisfactory feature since the perturbative- 
QCD corrections are known to be small. However a closer scrutiny makes it worrying since both 
the perturbative QCD corrections and the condensate predict a decrease which seems not to 
be seen. This leads us to start a very systematic study of the problem, with the following series 
of improvements on earlier works: 

• We reach an energy of 10 GeV by matching several lattice spacings so we are in a better 
position to eliminate lattice artifacts. 

• We make a systematic use of Ward identities relating the quark propagator and the vertex 
function via the constant Zy and study both quantities in parallel. 

• We make use of a very efficient way of eliminating hypercubic artifacts. A simpler version of 
it was elaborated while studying gluon propagators [Sj , . In this work we have encountered 
the necessity to improve significantly this method. 

This last point will be the main subject of this paper. Indeed, the raw data show a shape 
somewhat reminiscent of a half- fishbone (fig. ^) , utterly different from a smooth curve expected in 
the continuum. As we shall see the method elaborated in [7] proves not to be powerful enough. 
We therefore wish to attract attention on the generalisation of the above-mentioned method which 
we believe is strikingly efficient and should become a very useful tool for the lattice community. 

The remaining part of the work, i.e. the correction of S'0(4) symmetric artifacts and the 
resulting physics results will be presented in a later publication. 

In section 121 we will recall some theoretical premises, in section |21 we will indicate the lattice 
simulations which we have performed, in section|l]we will describe our method to eliminate lattice 
artifacts and compare it to earlier methods. 

2 Theoretical premisses 

We work in the Landau gauge. Let us first fix the notations that we will use. We will use all 
along the Euclidean metric. The continuum quark propagator is a 12 x 12 matrix S{p^). The 
inverse propagator is expanded 




(1) 



where a, b are the color indices. 
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Let us consider a colorless vector current q'y^q. The three point Green function is defined 

by 

G^^ip^q) = [ d^xd''ye'P-y-'^P+'^>^ < q{y)q{OhMm^) > ■ (2) 



In all this paper we will restrict ourselves to the case where the vector current carries a vanishing 
momentum transfer q^. The vertex function is then defined by 

r^{p,q = 0) = S-Hp)G^{p,q = 0)S-\p) . (3) 

In the following we will omit to write q/^ = and we will understand T^{p) as the bare vertex 
function computed on the lattice. The renormalised vertex function is ZvTf^{p). 
From Lorentz covariance and discrete symmetries 

r^(p) = Sa,b [giip^ht^ + ig2{p^)Pti + gsip'^)?)^^ + ig4{p'^)[yfi,^]] ■ (4) 

The Ward identity tells us that 

Zvr,{p) = -i^s-\p), (5) 



which from Q-© implies 



= Zygiip^), 2 -^Z^{p^) = Zygsip^), 
2^b{p'') = -Zvg2{p^), <?4(p') = 0. (6) 



For a conserved current, Zy = 1. We keep Zy since the local vector current on the lattice is 
not conserved; it will differ from 1 by lattice perturbative corrections which are a finite series in 
the "boosted" bare coupling constant, independent of p^. However, lattice artifacts do generate a 
sometimes significant p^ dependence of Zy at the level of raw lattice data, see for example fig. 01 
We will therefore sometimes use a p-dependent raw Zy written Zy{p^) and defined in eq. (fT7|) 

The renormalisation scheme that we use is the one called MOM', eq (26) in ref. [lUj . The 
bare propagator S{p) is multiplied by the renormalisation constant 



Sr{p) = Z^{fi)S{p), whence S^^{p)\p2^^2 = 5a,b (i;^ + m{p'^)) (7) 

'ipi 



This defines the quark field renormalisation. Due to the Ward identity, the factor Z^{^) ^ 



multiplies the bare vertex function gi, so that g^{p'^ = fj?) = Zy^. 

From the anomalous dimensions computed in ref. DJj we may express the perturbative 
running of Z^ for example as a function of the running o^joM (p) ■ in the continuum is a 

constant, Z^{p) and gi{p) have the same perturbative scale dependence. 
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' p=6.4 ; raw lattice data 

' result of the fit 



• f>=6.8 ; raw lattice data 

■ result of tile fit 



0.8 I ' • ' • ' • ' • ' • ' 1 0.8 I ' • ' • ' • ' • ' 1 

01234560 2 4 6 8 10 

p (GeV) p (GeV) 

Figure 1: Fit of the raw lattice data for clover- with formula for /3 = 6.4 and 6.8. For Q{p'^, 0, 0) 
we have taken a polynomial of degree eight in . The fit has been done independently for each /3. The 
raw data are represented by black circles while the fit is shown by red squares. The agreement is often 
such that the raw data and fit points are indistinguishable on the plot. 



3 Lattice calculations 

We have used improved Wilson quarks (often called clover) with the CSW coefficients computed 
in 100 gauge configurations have been computed at /? = 6.0,6.4,6.6,6.8 with volumes 24^, 
16^ and 8^. We have performed the calculation for five quark masses but in practice, for what is 
our concern in this paper, the quark mass dependence has non surprisingly proven to be negfigible 
and for simplicity we will only present the results for the lightest quark mass, about 50 MeV, i.e. 

K = 0.1346, 0.13538, 0.13515, 0.13489 for p = 6.0, 6.4, 6.6, 6.8 (8) 

It should also be mentioned that all the results presented refer to the 24^ lattices unless stated 
otherwise. 

We have also used overlap fermions jl2| I13j with about the same mass i.e. 

a mo = 0.03, 0.01667, 0.01. for /3 = 6.0, 6.4, 6.8 (9) 

with s = and volumes of only 16^ due to memory limitations. The bare mass mo and s are 
defined from 

Dover = {l + s + amo/2) + {l + s- amo/2) ^^"^"^^ + == (10) 

where D^{—{1 + s)) is the Wilson-Dirac operator with a (negative) mass term — 1 — s 

^w(-l -s)^ l^^iV^ + V* ) - JaV* -1-s, (11) 
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• P = 6.8 



0.5 1 1.5 2 2.5 3 3.5 4 2 3 4 5 6 7 

p (GeV) p (GeV) 



9 10 



2.8 



2.6 



Figure 2: Raw lattice data for overlap- Z.^ for /? = 6.0 and 6. 
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Figure 3: Raw lattice data for overlap- Zy for /3 = 6.0 and 6.8. 



The propagators S{x, 0) from the origin to point x have been computed and their Fourier 
transform 

S{p) = Y,e-'^-'=Six,0) (12) 

have been averaged among all configurations and all momenta within one orbit of the hy- 
percubic symmetry group of the lattice, exactly as for gluon Green functions in [Hj-|S]. In the 
case of overlap quarks the propagator is improved according to a standard procedure '13' which 
eliminates 0{a) discretization errors: 

SM = :r^^^ (13) 
1 — aniQ/Z 

From now on, the notation S{p) will represent the improved quark propagator in the case of 
overlap quarks and the standard one in the case of clover quarks. 
In both cases we fit the inverse quark propagator by 

S-Hp) = Kb Z^p^) (i^ + m(/)) (14) 
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p' extrapolation method 
• p'''' extrapolation method 



• p extrapolation method 
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Figure 4: comparison of the "pW extrapolation method, represented by black circles, with the " pl^"! 
extrapolation method " represented by "x" symbols. The left (right) plot shows the result for clover- 
at /3 = 6.4 {P = 6.8) 



according to eq. and where is defined in eq. (|18jl . The three point Green functions 
with vanishing momentum transfer are computed by averaging analogously over the thermalised 
configurations and the points in each orbit 



G^,{p, q = 0)=< 75 S{p)^ 75 7m S{p) > . 



(15) 



where the identity S{0,x) = 755"^ (x, 0)75 has been used. The vertex function is then computed 
according to eq. (jSJ and we choose for the lattice form factor gi : 



0) (7M-i'M^ 



(16) 



where the trace is understood over both color and Dirac indices. 

Finally, according to the Ward identity © we compute Zy simply from 



9i{p^ 



(17) 



where the p^-dependence of Zy coming from lattice artifacts has been explicitly written. 

In whole this paper we will use the values in the following table ^ for the lattice spacings, 
which follow the /? dependence found in ref. |15j . 
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6.0 


6.4 


6.6 


6.8 


a-' (GeV) 


1.966 


3.66 


4.744 


6.1 


a (fm) 


0.101 


0.055 


0.042 


0.033 



Table 1: Lattices spacings 
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Figure 5: Lattice data for overlap- after application of the " pP"] extrapolation method " for /? — 6.0 
and 6.8. 

4 Elimination of lattice hypercubic artifacts 

The question of eliminating lattice artifacts has been our main difficulty in the accurate study of 
the quark propagator. We became convinced that it was absolutely impossible to say anything 
sensible without an extremely careful elimination of artifacts. Here we mean mainly the ultraviolet 
artifacts, the infrared artifacts having never been really troublesome in this problem. 

We have elaborated a very powerful method to deal with hypercubic artifacts i.e. with those 
ultraviolet artifacts which come from the difference between the hypercubic geometry of the lattice 
and the fully hyper-spherically symmetric one of the continuum Euclidean space. The principle 
of this method is based on identifying the artifacts which are invariant for the symmetry of 
the hypercube, but not for the 50(4) symmetry of the continuum. 

Once these artifacts have been eliminated it is obvious, as we shall show, that other - 50(4)- 
invariant ultraviolet artifacts - are present. And these turn out to be even trickier to deal with, 
mainly because we did not fully understand their rationale. 

Therefore, we intend to restrict ourselves in this paper to a careful explanation of the hyper- 
cubic artifacts elimination method since we believe it represents a real progress and it can be 
useful for many other lattice calculations. The treatment of the 50(4)-invariant artifacts and of 
the physical results concerning the quark propagator will be given in a later publication. 

4.1 p[^"l extrapolation method 

Since we use hypercubic lattices our results are invariant for a discrete symmetry group, H4, a 
subgroup of the continuum Euclidean 50(4). This implies that lattice data for momenta which 
are not related by an H/j^ transformation but are by a 50(4) rotation will in principle differ. Of 
course this difference must vanish when a ^ but it must be considered among the discretisation 
effects, i.e. ultraviolet artifacts. For example, in perturbative lattice calculations one encounters 
the expressions 




(18) 
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Figure 6: Lattice data for overlap- after application of the " p[^"l extrapolation method " for (3 — 6.0 
and 6.8. 



Both are equal to up to lattice artifacts: 

p2 = V p2 = / _ l«2p[4] ^ . . . p2^p2_ 1 2^[4] ^ . . . ^ wherep^"! = V p^" . (19) 
r ^ o 

/i=l,4 A'=l,4 

All terms in the dots are proportional to a2np[2n-i-2] _ p2^ p2^ a?'^'^ are invariant under B./^ 
but only is under 50(4). For example the momenta 27r(l, 1, 1, 1) jL and 27r(2, 0,0,0) /L have 
the same but different and p^. In other words, if we call an orbit the set of momenta 

related by H/^ transformations, different orbits, corresponding to the same p^, will in general have 
different pt^' . The hypercubic artifacts can be detected by looking carefully for a given quantity 
at a given how it depends on the orbit. 

One method proposed with success for the gluon propagator [J] , |H] analyses a generic lattice 
measured quantity Q as a function (5(p^,pf^l). For a given value of p^, if enough different values 
of pt^l exist the quantity Q is fitted by 

where (5(p^,0) is free of hypercubic artifacts and where ^'^g^ '^^ |y=o is computed numerically for 
each from the slope of the lattice data for (5(p^,p''*^) as a function of p'^^. Of course, we could 
also consider pl^l, etc. but usually there are not enough different orbits for one given p^ to fit 
more than the pt^' correction. 

Let us call this method the "pW extrapolation method". For the gluon propagator this 
method has been shown [H] to lead to a resulting function G(p^,pW = 0) much smoother than the 
direct lattice results (5(p^,p'^'), even if the latter are restricted, as often done, to the "democratic" 
momenta, i.e. to those which have the smallest pW. This method could be applied with some 
success to the function clover-t^i (i.e. g\ obtained from improved Wilson quarks). 

But in general, as we shall see, when applied to clover-Z^ or to the quantities computed 
from overlap quarks, the "pt^' extrapolation" method fails. The signal of this failure is that the 
resulting function (5(p^,p[^^ = 0) still shows sizable oscillations typical of hypercubic artifacts. 



8 



1.03 



0.98 



0.93 - 



0.83 



1.05 



• result of p extrapolation method (P = 6.4) 
■ democratic orbits (|i = 6.4} 



^^i>. 3." 



'is 



3 

p(Gev) 



0.95 



0.9 - 



0.85 



• result of p extrapolation mettiod (P = 6.8) 
■ democratic orbits (|i = 6.8) 



••»........ 



p(Gcv) ^ 



10 



Figure 7: Red squares represent the "democratic" orbits, defined by p^^^/ip^Y < .5, for clover- Z,^ at 
P = 6.4,6.8. The black circles represent the result of the p'^"' extrapolation method. The latter exhibit a 
much smoother behavior. 



We then propose a "p[^"] extrapolation method" which allows to eliminate much more 
efficiently these hypercubic artifacts. The improvement goes in two directions: 

i) Instead of fitting the pt^l slope separately for each value of we try a global fit of the 
hypercubic artifacts over all values of 

ii) We chase hypercubic artifacts up to order a'^. 

In order to perform a global fit we start from the remark that in this paper we are dealing 
with dimensionless quantities, gi and Z^. It is thus natural to expect that hypercubic artifacts 
contribute via dimensionless quantities times a constant ^ . Next we assume that there is a regular 
continuum limit which implies that in the denominator we can have only physical quantities, 
namely ^ a function of p^. These two priors lead us to a Taylor expansion with terms of the type 

2k [2k+2n] \ 

k>0, n>0 k + n>l, m>0. (21) 

(^2)„ J 

This still leaves us with far too many terms to make sensible fits. It is reasonable to truncate 
this series in a and we choose to expand it up to a^. We will also truncate it to n < 1, and now 
comes an heuristic argument to justify this truncation. 

The lattice results are invariant and thus typically functions of p^, and p-p. Dimension- 
less quantities will depend on ratios of these quantities and on terms ^ like a^p^. An examination 
of these shows that they produce anisotropic terms corresponding to n = and n = 1 in the 
classification of eq. (|2H) . For example, a^p'^ produces from the expansion eq. a term cx a^p'^' 
while p ■ p/p^ produces terms oc a'^pW /p^ etc. This is also what is found in the free case jl3j . 

Hence we fit the data according to 

Q(/,aV^aV^---) = Q(/,0,0)+ci^ + C2 + C3 ^ + C4 a V^' (22) 

pZ \ P j P 

^We neglect a possible logarithmic dependence in p^. 
■^In all this discussion we consider the mass as negligible. 
^Terms like a'*(p ■ j3)p^ yield only 0[a^) anisotropic terms. 



9 



0.9 - 



0.4 



X ^=6.0 ; Vol=24 ; artifact-free result 
■ VqI=16' ; artifact-free result 

□ Vol=16' ; raw latticg data 



p(Gcv) 



1.05 



0.95 



0.9 



0.85 



3.4 



0.4 



^ ^=6.4 ; Vol=24 ; artifact-free resutt 

■ Vol=16' ; artifact-free result 
□ Vol=1 6^ ; raw lattice data 



!if, 



1.4 2.4 , 3.4 4.4 5.4 

p(Gcv) 



Figure 8: The white squares show the raw lattice resuhs for cfover-Z^ with a volume of 16^. The black 
squares are the same after the artifacts computed with a volume of 24^ have been subtracted and the 
crosses represent the artifact-free result computed with a volume of 24*. The agreement between black 
squares and crosses is striking except for the smallest momentum on 16*. The figure to the left (right) is 
for P = 6.0 {P = 6.4). 



A remark is needed about the a dependence of the coefficients q. Being dimensionless it 
is expected in perturbation theory that these coefficients depend only logarithmicahy on a and 
taking them as constants would seem reasonable. This conjecture does not work as shown in fig. 
ini which is the sign of non-perturbative ©(aAqcD) contributions. Still this figure shows a rather 
convincing linear dependence of the q's for which tells that a good global fit can be performed 
by expanding the coefficients: Ci{a) = + acj. For gi the a dependence is not linear while the 
hypercubic artifacts are one order of magnitude smaller. 



The functional form used for Q{p'^, 0,0) does not influence significantly the resulting artifact 
coefficients. We can even avoid using any assumption about this functional form by taking all 
the values for 0, 0) as parameters which can be fitted ^. 

This improved correction of hypercubic artifacts turned out to be particularly necessary for 
Z.^. In fig.^the very strong hypercubic artifacts produce an impressive branched structure with a 
kind of periodicity. In fig. 0]we show the effect of both the use of eqs. ()2U() and ()22() . It turns out 
to that "pW extrapolation method", eq. H2U() . makes the branches disappear and the curve 
look much smoother. However it still contains some oscillations reminiscent of the hypercubic 
artifacts. The "pP"l extrapolation method", eq. 1)221) . brings in a further dramatic smoothing. 

The same is true for overlap-computed quantities. In figs. ElandOlthe raw lattice data for Z^ 
and Zv exhibit dramatically the "half-fishbone" structure which is a symptom of strong hyper- 
cubic artifacts. In figs. El and El the same data are shown after applying the pl^'^l extrapolation 
method. Clearly the curves are now perfectly smooth. We will return later to the fact that Zy 
is not a constant. 

Altogether we would like to stress the following hierarchy: First, the hypercubic artifacts are 
one order of magnitude larger for overlap quarks than for clover ones, see fig. El Second, for both 

*We have enough data for that. 
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Figure 9: CoefScients of the hypercubic artifacts for from overlap quarks (left) and clover (right) as a 
function of a for /3 = 6.0,6.4,6.8 and 6.6 in the clover case. The squares corresponds to ci, the coefficient 
of a'^pW /p^, the triangles to C2, the coefficient of (a^pW /p^)^, the stars to C3, the coefficient of (a^'p'^l /p^)^ 
and the circles to C4, the coefficient of a^pW. It suggests a linear dependence on a, especially for the 
overlap quarks. The overlap artifacts are larger than the clover ones by one order of magnitude. 



types of quarks the hypercubic artifacts for are one order of magnitude larger than those for 
91- 

4.2 Comparison with the "democratic" method 

The hypercubic artifacts, sometimes called "anisotropy artifacts" have been a long standing 
problem in lattice calculations. Studying the gluon propagator the authors of ref. [Hj where 
aware that the problem was related to the fact that for a given momentum these artifacts were 
minimised when the components were as small as possible, i.e. such that the components are not 
too hierarchical, and that the ideal situation was the diagonal p oc (1, 1, 1, 1), whence the name 
commonly used of a "democratic" repartition of the momentum in all directions. Therefore they 
have proposed a selection keeping only the orbits having a point within a cylinder around the 
diagonal. Several other criteria have been used. 

In this subsection we want to compare this method of eliminating the non-democratic points 
to the " a^pl-^"! extrapolation method", eq. ()22|). If we try to select, Jl], the orbits which are in 
a cylinder around the diagonal with a radius 2it/L, we are left with only 11 orbits among 69. 

In order to have a less restrictive criterion and to make the bridge with the method used here 
we will use the pt'^'^l's defined in eq. (|19|) . In our language, democracy can be translated as a 
small enough ratio /(p^)-^. Momenta proportional to (1, 1, 1, 1) and (1, 0, 0, 0) have ratios 1/4 
(minimum) and 1 (maximum) respectively. In fig. [3 we plot for the result of the following fit. 
We take the "democratic" orbits defined by ^[^^/(p^)^ < .5. This leaves 40 orbits out of 69 for 
every p. Fig. [7| clearly shows oscillations demonstrating that the hypercubic artifacts have not 
been totally eliminated. For this reason and also because of the loss of information due to the 
rejection of "undemocratic" points, we did not use this method. 
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4.3 Finite volume artifacts 

We did not see any sizable finite volume effect in the case of clover quarks. To illustrate this claim 
we have performed the following exercise illustrated in fig. |H1 We have subtracted from the raw 
lattice results clover- Z^, computed with a volume of 16^, the artifacts with the coefficients ci, .., C4 
fitted on a volume 24^ , namely the results of eq. H22() and compared the result to the artifact-free 
function (5(p^,0,0) computed with 24^. The agreement as shown in fig. El is impressive except 
for the smallest momentum on 16^. We have also checked on several examples that the inclusion 
in the fits of finite volume artifacts of the type l/(L^p^) did not produce any significant change 
in the results. 

4.4 Zy and other discretization artifacts 

Figure El presents the result of extracting the hypercubic artifacts from the raw lattice data for 
Zv computed according to eq. p7)) . It shows that we are not through with artifacts. Indeed, 
as we have already mentioned, the artifact-free Zy must not depend on p^. It is expected to 
depend on the bare coupling constant i.e. on /3 but not on the momentum. The figures El show 
smooth curves which confirms the efficient elimination of hypercubic artifacts, but it also shows 
a residual significant dependence on up to 50 % variation in the case of overlap quarks at 
P = 6.0. This dependence is necessarily due, either to additional finite lattice spacing artifacts 
which are not of the hypercubic type but are S'0(4)-invariant, or to finite volume artifacts. It is 
however noticeable that at /? = 6.8 Zy is really flat except for the two flrst points. 

These artifacts do not look simple since Zy increases at small p'^. We have discarded flnite 
volume effects in the preceding section for clover quarks but we could not, by lack of computing 
resources, perform the same check for overlap quarks. In particular the strong dependence of 
Zy at small momentum seen in fig. (HI is evocative of oc l/(L^p^) finite volume effects. But these 
type of effects would produce exactly the same shape at both /3's. The difference between the 
two plots in fig. shows that the case is more subtle and small momentum ultraviolet artifacts 
can also be present. This clearly needs a careful study and some theoretical understanding which 
will be developed elsewhere. 

Finally it is useful to notice that the p'^ dependence of Zy is one order of magnitude larger 
for overlap quarks than for clover ones. 

5 Discussion and conclusions 

We have computed the quark field renormalisation constant Z^ and the vector current form factor 
gi{0) both with improved Wilson quarks (clover) and with overlap quarks. The quark propagator 
is very strongly affected by lattice artifacts ^. This is already very annoying in the case of clover 
quarks but is even one order of magnitude larger for overlap quarks. See for example figures ^ 
and El Hypercubic artifacts mainly affect Z^, while the vertex factor gi is less affected by one 
order of magnitude. This forwards an invitation to use intensively Ward identities and vertex 
functions simultaneously to the propagator. 

In order to eliminate the latter hypercubic artifacts we have improved the method presented 
in refs. [71 into what we call the "pl-^"! extrapolation method". It is based on a systematic 
expansion over the invariants of the hypercubic group H/^^ which are not invariants of the 5*0 (4) 

^Notice that the size of lattice artifacts might depend strongly on the parameter s defined in ea.po^ and which 
we have taken to vanish in this paper. 
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symmetry group of the Euclidean continuum and on a systematic use of dimensional arguments 
to guess the dependence of the artifacts. We have shown that this method totally eliminates 
the dramatic disorder of the raw data, see figures which exhibit a shape vaguely reminiscent 
of the half of a fishbone. 

After applying the "p[^"l extrapolation method" we get results which are perfectly smooth, 
figs.|l]|ni In particular the artifact fig. [5] for overlap quarks should be compared to the raw lattice 
data fig. 121 In fig. ^ there is also a very good agreement between the clover raw data (black 
circles) and the red squares computed from eq. ()22() . We have shown that this method is much 
more efficient than the popular "democratic" one, and it is also more systematic and allows the 
use of all the lattice data which improves the statistics. We have decided to center this paper on 
this method because, although the point might look technical, we believe that it can be of great 
help for the lattice community when artifacts are large. 

The next task is to find out an equally efficient method to eliminate the isotropic arti- 
facts. Their presence is obvious from fig. El which shows a strong dependence of Zy after 
the anisotropic artifacts have been eliminated by the p[^"] extrapolation while Zy should be a 
constant ^ . Contrarily to Zy , Z^ = Zy gi is expected to depend on p'^ as an effect of the per- 
turbative QCD running and of the nonperturbative < > condensate. Interesting physics can 
be learned from this dependence provided we manage to fully control the isotropic artifacts. The 
constancy of Zy, being a strong constraint, will be a significant check that this control has been 
achieved. 
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